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Abstract 

We discuss the effects of the solvent composition on the helix-coil 
transition of a polypeptide chain. We use a simple model to demon- 
strate that improving the hydrogen-bonding ability of the solvent can 
make the transition less cooperative, without affecting the transition 
temperature. This effect is very different from other solvent effects 
which primarily influence the melting transition rather than the co- 
operativity. 
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Proteins and nucleic acids (DNA and RNA) are biopolymers which are 
involved in almost every biological process. These molecules are character- 
ized by special hierarchical structures Each chain consists of a unique 
sequence of bases (DNA) or amino acids (proteins), called the primary struc- 
ture, which determines its specific biological function. Short-range interac- 
tions between nearby monomers along the chain can lead to the formation 
of a three-dimensional structure, usually a helix or a planar sheet, which 
is called the secondary structure of the chain. The secondary structure is in 
large measure determined by the primary structure. However, it also depends 
on the environmental conditions that the polymer experiences like the tem- 
perature or the composition of the solvent in which the chain is embedded. It 
is well known, for instance, that upon increasing the temperature, the helix 
structure "melts" into a random coil structure ||. Another example is the 
folding-unfolding transition taking place when proteins are stretched by ex- 
ternal forces 0]. In this paper we study the effects of the solvent composition 
on the helix-coil transition with the emphasis on a particular effect occurring 
when hydrogen-bonding agents are added to the solvent. We present a simple 
model which quantifies this effect, and compare it to other solvent effects, 
as well as to the effect that temperature or externally applied force have on 
the transition. While we will mainly deal with helix-coil transition occurring 
in polypeptides, we believe that the model may also apply to the transition 
observed in DNA experiments, just as the temperature-induced transition in 
both systems is studied by means of the same model. 

Polypeptides are chains of amino acids monomers linked to each other 
by covalent peptide bonds. Each monomer can exist either in a a-helical 
(h) or a coil (c) state. In a helical domain hydrogen bonding occurs along 
the polypeptide backbone between monomers separated by approximately 
4 units. These hydrogen bonds make the h state energetically more favor- 
able than the c state. The c state is entropically favorable because of the 
rotational degrees of freedom of the amino acids. The helical monomers at 
the boundaries between h and c domains represent a special case: These 
monomers lose their configurational entropy with no energy gain due to the 
formation of hydrogen bonds. The thermodynamic behavior of the chain 
can be described using the Zimm-Bragg model [|J, in which the free energy 
associated with a particular configuration (a sequence of h and c monomers) 
of the chain is given by 

F conf = NhAF + 2N d AW. (1) 
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In the above equation A^ denotes the number of monomers in the h state, A^ 
is the number of helical domains, AF is the excess free energy per monomer 
in the helical state, and AW is the cost in free energy to create a bound- 
ary separating h and c domains. The partition function is readily obtained 
using the transfer matrix approach which is also applied to solve the one- 
dimensional Ising model ||. In the Zimm-Bragg model, the transition matrix 
is given by 

M-( M — M ^)-( 1 (2 ) 



where 

s = exp(-AF/A;T), (3) 

and 

a = ex.p(-2AW/kT), (4) 

are the so-called Zimm-Bragg parameters, k is the Boltzmann constant, and 
T is the temperature. At low temperatures AF < (s > 1) since the energy 
gain in having an amino acid at a A state overcomes the entropy loss compared 
to the c state. Therefore in the low T regime most of the monomers are in 
the h state. At high temperatures s < 1, and the chain is mostly in the c 
state. The helix-coil melting transition occurs when s ~ 1. In contrast to 
the significant variation of s with the temperature, the parameter a (also 
called the cooperativity parameter) has only weak dependence on T. For 
biopolymers a is typically very small in the range 10~ 3 — 10~ 4 . The smallness 
of cr reflects the large free energy penalty in creating a helix-coil interface, 
and is the reason for the sharpness of the helix-coil transition (which, due to 
the one-dimensional nature of the system, is not a true phase transition). 

To understand the helix-coil transition in a more quantitative manner, 
we write the free energy of a chain consisting of A^ (N ^> 1) monomers: 

F ~ —NkThaxi, (5) 

where x\ is the largest eigenvalue of the transfer matrix M [Eq.(^]. The 
eigenvalues of M are the roots of the quadratic equation 

(1 — x)(s — x) — sa = 0. (6) 

The largest root of this equation is equal to 

' (7) 
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Figure 1: The fraction of monomers in the h state, ph = (Nh)/N, as a 
function of the Zimm-Bragg parameter s, for a = 10~ 3 (dashed line) and 
a = 1 (solid line). The curves intersect at s — 1, where ph = 1/2. 



The average number (Nh) of monomers in the h state is derived by differen- 
tiating the free energy [Eqs.@ and ([?[)] with respect to AF = —kTlns. We 
obtain 
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The dependence of ph = (Nh)/N, the fraction of monomers in the h state, on 
s is plotted in Fig. [IL for o = and a = 1 (the non-cooperative case). We 
observe that ph ~ for s Cl, while ~ 1 for s ^> 1 . The crossover between 
these two regimes occur at s ~ 1 (at s = 1, the fraction of helical monomers 
is exactly 1/2). The transition becomes sharper with decreasing a. For small 
values of a the width of the transition regime scales as As ~ o 1 ! 2 . 

When dealing with the possible influence of the solvent composition on 
the transition we should account for two types of effects. Effects of the first 
type decrease the relative stability of the h state over the c state. Effects of 
the second type reduce the helix-coil interfacial free energy. 

Within the first group, we include changes in the pH or the ionic strength 
of the solution that destabilize the h state 0. One can incorporate these 
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changes in the Zimm-Bragg model by a proper redefinition of the parameter 
s, which now becomes a function of both the temperature and the solution 
conditions. As a result of the change in the solvent conditions, the transition 
(melting) temperature shifts, and it is now defined by s*(T m ) = 1 with s* 
being the solvent-dependent parameter. 

It is interesting to note that the force-induced transition, observed in 
recent stretching experiments of DNA Jj]], has been also explained by the 
Zimm-Bragg theory with a rescaled (in this force- dependent) param- 

eter s ||. The effect of externally applied force on the transition is usually 
studied in the fixed-force (/) ensemble. On a mean-field level (which turns 
out to provide a rather good description of the experimentally measured 
elastic behavior), the elastic free energy is given by 

F d (f) = N h Ettf) + (N- N h )E?(f), (8) 

where E^(f) and Ef(f) denote, respectively, the elastic free energy per 
monomer in the pure h and c states (i.e., when all the monomers are in the 
same state). The elastic behavior of the pure phases is usually described by 
one of the generic polymer models like the freely-jointed or worm-like chain 
models |§ . Within these models, the elastic energy per monomer is given by 

E e \f) = -^J Q f R(f')df>, 

where R{f) is the (mean) end-to-end separation of the chain. The transfer 
matrix corresponding to the sum of configurational (|I|) and elastic (|5p free 
energies is 

and by comparison with the Zimm-Bragg transfer matrix (|2|) we readily con- 
elude that s*{f) = se K el (/)-^(/)]/*T 

Let us consider the change in the behavior of polypeptide chains that 
occurs when the hydrogen-bonding ability of the solvent is improved. One 
may consider this as another example of a solvent effect of the first type. In 
alcohol environment, for instance, the hydroxide groups tend to form hydro- 
gen bonds with the coil monomers, thus increasing the relative stability of 
the latter compared to the helical monomers. This effect is enhanced if the 
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hydrogen-bonding ability of the solvent with the c monomers is improved, 
e.g., if water is added to the solution |10|. However, one may consider the 
following different scenario which serves as an example of a solvent effect of 
the second type. Suppose the polypeptide chain is dissolved in an aqueous 
solution, and most of the coil monomers are involved in hydrogen bonds with 
the water molecules. Let us assume now that the hydrogen-bonding ability 
of the solvent is improved by adding molecules which prefer to bind to the 
terminal monomers of the helical domains. The high affinity of the new 
hydrogen-bonding agents (which we shall call "impurities") for the terminal 
h monomers leads to the formation of new hydrogen bonds at the interfaces 
between h and c domains. To introduce the effect of these additional hydro- 
gen bonds, we need to add a term to the configurational free energy [Eq.(|TJ)] 

F conf = NhAp + 2N d AW + /iN h (9) 

where /i is the chemical potential per impurity molecule attached to the chain 
at the interface between h and c domains, while N[ (iVj < 2Nd) denotes the 
number of such impurities [11|. It is useful to write 

/i = AW - AW - kT In (p, (10) 

where AW is the free energy per "impure" interface with an additional hy- 
drogen bond (replacing AW, the free energy of a "pure" interface), while 
ip = N£ 3 /V is the concentration of the free impurity molecules in the solvent 
(with £, the "thermal wavelength"). The last term on the r.h.s of Eq. fllOD 
represents the loss of mixing entropy of the impurity molecules connected 
to the chain compared to the free impurities (assuming that tp is sufficiently 
low). All the other contributions to the free energy of the impure interfaces, 
for instance the binding energy between the impurities and the chain, are 
included in AW. 

In order to find the free energy corresponding to Eq. (Q) we define two new 
states h* and c* representing, respectively, the terminal ("C-end") monomer 
of a helical or a coil domain followed by an impure interface. We can then 
apply the transfer matrix approach, where in the present case the transfer 
matrix is the 4x4 matrix 
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with 



—2 [AW — kT\nip) jkT = exp (-2AW/kT) <p 2 . 



a = exp 

The eigenvalues of the transfer matrix ([TT|) are the roots of the polynomial 
equation 

"{l-x)(s-x) -s(a + 2Vaa + a) = 0. (12) 
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Although this is a fourth order equation, the largest root (which is relevant 
to the thermodynamic behavior) comes from the quadratic equation defined 
in brackets in Eq. (|T2"D 



(1 - x)(s - x) - s(a + 2\faa + a) = 0. 

This equation is similar to the quadratic equation of the Zimm-Bragg model 
[Eq.(^|)], with the cooperativity parameter a replaced by 

\fo* = yfo + \fd = \fo + exp {—AWjkT\ (p. (13) 



Our model, therefore, yields a description of the helix-coil transition similar 
to the Zimm-Bragg model, but with a larger cooperativity parameter. We 
find that a* > a even if the chemical potential /i in Eqs.([|) and ( |I0D is 
positive, i.e., if the binding of the impurities to the chain increases the free 
energy of the system. This result is due to the simple fact that some of the 
helix-coil interfaces will have impurities even when [i > and, consequently, 
the free energy cost per interface will decrease (compared to the case when 
we had no impurities at all, <p = 0). In the "neutral" case \x = 0, for instance, 
half of the hydrogen-bonding impurities will be attached to the chain, and 
from Eqs.([|), (0), and (0) we find that a* = 4a. We thus conclude that 
the helix-coil transition becomes broader in the presence of the hydrogen- 
bonding impurities, but quite unusually, the melting temperature (which is 
associated with the parameter s) is not predicted to change. This effect is, 
therefore, markedly different from the other ("first type") solvent effects, as 
well as from the effect of externally applied force on the transition. 

To summarize, we propose a simple description for the complicated de- 
pendence of the helix-coil transition on the solvent character. We suggest 
that solvent effects can be characterized by their impact on the Zimm-Bragg 
parameters. Accordingly, we can broadly classify them into two types: Sol- 
vent effects of the first type are associated with a shift in the position of the 
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helix-coil transition, and the rescaling of the parameter s. Solvent effects of 
the second type reduce the cooperativity of the transition, i.e., increase the 
parameter a. The amount by which o grows depends on the binding free 
energy of the impurities, and is proportional to their concentration in the so- 
lution. This, however, applies only to low concentrations. At larger densities 
it is necessary to include higher virial coefficients in the mixing-entropy term 
[Eq JlOD l to describe the interactions between the impurities. 
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